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Abstract 

We analyze various models for the top quark forward-backward asymmetry (A‘,,) at the Teva- 
tron, using the latest CDF measurements on different Aj,,s and the total cross section. The 
axigluon model in Ref. has difficulties in explaining the large rapidity dependent asymmetry 
and mass dependent asymmetry simultaneously and the parameter space relevant to A‘,, is ruled 
out by the latest dijet search at ATLAS. In contrast to Ref. , we demonstrate that the large 
parameter space in this model with a U(1)q flavor symemtry is not ruled out by flavor physics. 
The t-channel flavor-violating Z’, W’ and diquark models all have parameter regions that satisfy 
different Arg measurements within 1 ø. However, the heavy Z’ model which can be marginally 
consistent with the total cross section is severely constrained by the Tevatron direct search of same- 
sign top quark pair. The diquark model suffers from too large total cross section and is difficult 
to fit the tt invariant mass distribution. The electroweak precision constraints on the W’ model 
based on Z’-Z mixings is estimated and the result is rather weak (mz > 450 GeV). Therefore, the 
heavy W’ model seems to give the best fit for all the measurements. The W’ model predicts the 
tt + j signal from tW” production and is 10%-50% of SM tt at the 7 TeV LHC. Such t+ j resonance 


can serve as the direct test of the W’ model. 
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The prompt decay of top quark before hadronization provides opportunity to explore 
its various properties like charge, mass and spin. Given its large mass, the scale of top 
quark pair production is greater than 2m, where the perturbative QCD plays important 
role. Therefore, the top quark pair production at hadron colliders can serve as handle of 
precision test of the standard model (SM) gauge interaction, both weak interaction in its 
decay and the perturbative QCD theory of strong interaction in its production. 

From the structure of the SM, top quark is special. As a colored particle, it is the 
heaviest known particle which is copiously produced at the hadron collider. Since the top 
quark acquire its large mass through the electroweak symmetry breaking (EWSB), any of its 
properties deviated from SM would be an important signals for new physics and potentially 
indicate the origin of EWSB, which makes searching new physics in the top quark sector 
extremely interesting at both Tevatron and LHC. 

One important measurement for top quark in top quark pair production is the top forward 
backward asymmetry, which is equivalent to charge asymmetry under CP transformation 
. For the SM production, it involves the high precision calculation of QCD. At O(a), the 
bremsstrahlung amplitudes 99 — QQg carry an odd power of color charge hence have an 
odd charge conjugation parity in the interference terms among initial states radiation and 
final state radiation diagrams. There is also interference between the box diagram of O(a‘) 
with the LO diagram that contributes to the charge asymmetry. 

CDF collaboration has recently updated the measurements on the total forward-backward 
asymmetry in top quark pair production with the semi-leptonic tt data with integrated 
luminosity of 5.3 fb™t |}. The observed total asymmetry measured in the lab frame and the 


tt rest frame are 

Atp = 0.150 + 0.050(stat) + 0.024 (syst) (pp rest frame) 

Atp = 0.158 + 0.072(stat) + 0.017(syst) (tt rest frame) (1) 
which corresponds to the SM prediction based on the NLO simulation, Monte Carlo for 
FeMtobarn processes (MCFM), 0.038 + 0.006 (in lab) and 0.058 + 0.009 in tt rest frame 


respectively |3]. These measurements have improved the previous results based on 3.2 fb=! 


of APP,(cos9) = 0.19 + 0.069 and A%,(An) = 0.24 0.014. [3] 2 


1 The top quark forward-backward asymmetry has also been measured in the di-lepton channel as Arp = 


0.42 + 0.15(stat) + 0.05(syst) in the tt rest frame with 5.1 fb~+ data [2]. 
? Note that the recent DO measurement At p = (8 + 4(stat) + 1(syst))% is based on top-pair events that 
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More importantly, with the enlarged data sample, CDF collaboration has also released 
two distributional measurements. The most interesting result is the mass dependent forward 


backward asymmetry. The mass dependent forward backward asymmetry in the tt rest frame 
Ali „(Mir > 450 GeV) = 0.475 + 0.112 (2) 


in comparison to the QCD correction prediction 0.088 + 0.013. This 3.5 ø deviation may 
be a strong indication for physics beyond the SM. The second measurement is the rapidity 


dependent asymmetry, which is frame independent, as 
Ars(| Ay |> 1.0) = 0.611 + 0.256 (3) 
Ara(| Ay |< 1.0) = 0.026 + 0.104 + 0.056 (4) 


in comparison to the MCFM SM prediction as Arg(| Ay |> 1.0) = 0.123 + 0.018 and 
Arg(| Ay |< 1.0) = 0.039 + 0.006. 

The ratio of the parton level asymmetries in the two different frames, which differ by 

longitudinal boost, is 

App 

Gr = 0.95 £0.41 (5) 
with the error corrected for the expected correlation across frames in the NLO QCD assump- 
tion. Even though the uncertainty is still large, this close to 1 central value implies that the 
top events which contribute to the asymmetry mostly lie in the forward-backward direction 
so the asymmetries are less dependent of the longitudinal boosts along the beam direction. 
This feature is also shown in the An dependent asymmetry A (| Ay |> 1.0) = 0.611 which 
shows that the asymmetric events are mostly due to events with larger rapidity difference 
Im — m |: 

On the other hand, the measurement of tt cross section o, updated by the 4.6 fb~' CDF 
result (with m, = 172.5 GeV), is a7? = 7.50 + 0.31(stat) + 0.34(syst) + 0.15(Z theory) pb 
which is in very good agreement with SM theory prediction of ot? = 7.5105 pb at NNLO Ël. 
Therefore, in order for new physics to generate large asymmetry without changing the total 


production cross section, the new physics contribution must interfere with the leading SM 


satisfy the experimental acceptance, which is uncorrected for effects from reconstruction or selection and 


can not be used to compare with the CDF results 4. 
3 The latest NNLL calculation shows the o,¢(m; = 173.1 GeV) = 6.30+0.19+0:31 pb [1] which is significantly 


lower than the experimental results. However, we still use the old SM predictions since we do not know 


clmi = 172.5 GeV) for the latest results. 


production of wu, dd + té as color octet exchange in s-channel. For instance, in order for the 
s-channel massive Z’ to explain the asymmetry, there is no interference between s-channel 
color singlet exchange ut, dd Z, tf and QCD uu,dd + tt. The asymmetry events due to 
Z' will significantly enhance the total cross section g at the same time and this causes a 
strong tension between fitting of Arg and og. This requirement implies that there are only 


two categories of candidate models to solve this anomaly. 


e First category of models contain s-channel color octet vector boson but with parity 


violation at both q — ¢— G and t — t — G vertices 


e Second category correspond to the t-channel exchange of light gauge boson of maximal 
flavor violation that couples initial state u,d quark to the third generation 1 
Whe 


The large asymmetry can be generated via Rutherford singularity behavior 


Both categories of models have their realizations in the beyond SM models. Given the 
updated measurements, especially the new distributional measurements, we discuss the cur- 
rent status of various models. In addition, the models may have other implications that 
have been or will be constrained by some direct or indirect experiments. One realization 
of the first category models is the non-universal axigluon model proposed in [öl and it may 
receive constrain from low energy neutral meson mixings |8]. However, we show that the 
flavor bound can be easily evaded by putting a horizontal flavor symmetry U(1)g. In the 
Wr models hd, since the Wr is charged under SM U(1)em, the neutral component wè 
would inevitably mix with W? and some extra U(1)x which induce a Z-Z’ mixing. The new 
ATLAS Dijets search and the Tevatron same sign dileptons [18] would severely constrain 
the s-channel axigluon models and the t-channel heavy Z’ models. We also study the direct 
prediction at the Large Hadron Collider(LHC) using the 1 o fitting of all three asymmetry 
measurements A% „(Mir > 450 GeV), A%,(|Ay| > 1.0), A%,(total) with the right total cross 
section. 

The paper is organized as follows. In Section [I| we presented the 1 o fitting of all 
three asymmetry measurements for s-channel color octet model (Section [LA], t-channel 
Z' model (Section [[B1), W’ model (Section [[B 2), diquark model (Section [[B3) and the 
corresponding consequences. In Section we calculate the production rates for the new 
particles in various different models at the Tevatron which give the bounds for those models 


and the LHC signals. In Section [LI] we consider some indirect bounds for the axigluon 


4 


from flavor physics (Section IA) and W” model from electroweak precision test (EWPT) 
(Section MIB). Section [V] contains our conclusions. 


I. UPDATED STATUS OF THE MODELS 


In this section, we discuss the updated status of the models based on the latest mea- 
surements, especially the new distributional measurements. The simulation in the following 
discussion is at parton level and leading order. The asymmetry observables are defined at 
parton level without taking into account possible reconstruction efficiency. The SM contribu- 
tion to the asymmetries from MCFM simulation have been subtracted to the corresponding 
measured values. The total cross section is obtained by multiplying a QCD k-factor. Since 
the latest experimental value is based on m = 172.5 GeV, for better comparison, we em- 
ploy the theory calculation at NNLO for m; = 172.5 GeV and the k-factor is 1.3. Last, the 
differential cross section of tt invariant mass is not included as requirement in the scan since 
QCD correction and cut efficiency may significantly modify the shape of differential 
distribution da/dy,,;. 

In the following discussion, we are mostly interested in the region where the three asym- 


metry measurements can be explained within 1 ø. 


A. s-channel color octet 


The interference term between the color octet V — A gauge boson Gi contribution and 
gluon contribution in qq annihilation pickup a term as 
2958(8 — Ma) 
(8 — M@)? + Male 


where gs is the strong coupling, g“ is the axial component of the coupling between Gi, and 


[+2 9% g% 3 cos 6] (6) 


light quarks q and g‘, is that of the top quark. If the interference contribution is positive 
asymmetry, it requires that the axial coupling g4g', < 0 is inevitable |. 
The s-channel models can be realized in various context. The first realization is the 


axigluon models where SU (3). color gauge symmetry is only a remnant of SU(3)Lx SU(3)r 


4 This non-universal gauge interaction potentially cause the violation of GIM mechanism thus may be 
constrained from flavor changing neutral current (FCNC) processes such as neutral meson mixings and 


we discuss its implications in the next section. 


broken by a bi-triplet scalar and another color octet with axial coupling become massive. 
However, to achieve the g4g‘, < 0 requirement, the axigluon model has to be non-universal 
and one example is the 4-generation model proposed in N Another realization is the models 
of extra dimension theory where massive color octet Kaluza-Klein (KK) gluon couple to the 
SM quarks in the chiral form as a result of fermion profiles . The large m, naturally 
implies that the top quark and light quarks couple to KK gluon in different way. 

One interesting feature that was discussed in the axigluon model |] is the mass dependent 
asymmetry. Due to the opposite contribution to asymmetry between the interference term 
and new physics squared term, the asymmetry is positive when the centre-of-mass energy 
is at intermediate energy but when it is close to the threshold, the asymmetry may become 
negative. This bending-over in correlation between asymmetry A’ and the centre-of-mass 
energy M had been shown in the latest CDF measurements, in both the measurement with 
finite bin sizes of Mz and the measurement with below/above Miz edge. 

We use the axigluon model as one example to illustrate the feature of s-channel models 


in comparison with the updated measurement. Figure [I| shows the summary of best fit 
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FIG. 1. 1 o parameter region for constraints from the total asymmetry in the lab frame and the tt 


production rate; 1.5 o parameter region for the Arg in Ma > 450 and | An |> 1 and o(pp — tt). 


parameter regions for total asymmetry, mass dependent asymmetry, rapidity dependent 


asymmetry, the total cross section and the last bin of do/dM;r measurements. Since the 
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total asymmetry has been reduced from the previous fitting in I the 1 o region with total 
asymmetry is enlarged as shown in Fig[i] However, there is no 1 ø region for the mass 
dependent asymmetry of Mq > 450 GeV or the rapidity dependent asymmetry of | An |> 1. 
Figure [I] shows the 1.5 o parameter space for Arg(| An |> 1), Arg(M;r > 450 GeV) as well 
as the total tt production rate o(pp — tt). It is clearly shown that the axigluon model 
does not consistently generate the large asymmetries in the events of Mẹ > 450 GeV and 


An j> 104 


B. t-channel 


As we argued, the ratio of APP /A' close to one may imply that the top events are mostly 
in the forward-backward direction so the asymmetries are less dependent of the longitudinal 
boosts along the beam direction. Since the t-channel models naturally predict large number 
of events in the forward-backward region, the close to one ratio of A??/A" is a basic feature 
of t-channel models. 

If the asymmetry is due to new physics in t-channel physics, the interference contribution 


between new physics and SM QCD is proportional to 


ISIN 
Or (u? + 8m? +...), (7) 
Sty 
where t; = —35(1 — 6cos@) and 1/t expansion naturally picks up a cos@. The t-channel 


physics naturally generates a large asymmetry in the tt system. In addition, the maximal 
asymmetry is generate at the Rutherford singularity where 9 = 0 which corresponds to very 
high centre-of-mass energy. One would then expect the positive correlation between A‘, 
and Ma. 

The t-channel Z’ model in ld proposed a color singlet neutral gauge boson with max- 
imal flavor violation between first and third generations and the new contribution interferes 
with the SM wu — tt. Similar to the Z’ model, instead of neutral current exchange in 
t-channel, there is also a proposal using charge current exchange in t-channel as flavor vio- 


lation W’. The interference effect is reduced since it’s only the dd initial state hd . Such 


5 One can use the general color octet vector boson with V — A interaction to fit the two distributional 
asymmetries and total asymmetry [3]. The general results together with the most recent bounds in the 


dijet channel from ATLAS will be presented elsewhere. 


flavor violation gauge interactions may be realized in horizontal gauge interaction models 


for neutral current or generalized left-right model 


for charged current. 
A Higgs-like scalar with maximal flavor violation|22 would generate a large negative 
asymmetry due to the helicity-flip in the Yukawa coupling. The spin conservation in the 
0 = 0 direction requires the top quark to move backward. To resolve this, the fermion- 
number violating diquark scalars with maximal flavor violation was proposed hul Ld m 
Diquark scalar can be 3 & 3 = 6 3 under SU(3). and has fermion-number violating 
coupling as tuo or td@. Such diquark scalars with flavor violation can also be realized 
in various BSM contexts, partial unification models or supersymmetry. For instance, R- 
parity violation supersymmetric standard model which contains the baryon number violating 
coupling, Eugyusdzd} mt the down type squark d; can mediated u-channel dd — tt that 
interferes with the QCD dd — tt. 

All the three proposals can in principle predict large positive asymmetry in tt production. 


In the following paragraphs, we examine the numerics to see whether the models can explain 


the three asymmetry measurements and the total cross section at the same time. 
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We first examine the first proposed t-channel model, Z’ l, To minimize the constraints 
from low energy, the authors proposed a right-handed coupled Z’ with large coupling between 
u and t. The parameter region for 1 o fitting of all three asymmetry measurements as well 
as the total cross section for light Z’ mostly below t threshold is presented in Figure P} Due 
to large destructive interference, the total cross section is always smaller than the measured 
value. This result is also shown in the NLO calculation of Z’ model Ri The best fit points 
for heavy Z’ by requiring 1 o fitting for all the three asymmetry measurements are listed 
in the Table [I| The corresponding tt cross section are also below the 1 ø total cross section 
and the best points are towards heavy masses of O(700 GeV). 

One more complication which has been discussed in is that the events in the 
t-channel exchange tend to be in high energy region which significantly increase the tail 
of da/dMz, especially the last bin (800 GeV-1.4 TeV) in da/dM,. The QCD correction 
may change the shape and lower the contribution in high energy mi In addition, the t- 


channel kinematics implies that the top quark events at high energy are mostly in the larger 


0.7 ——— 


0.6) 


0.5? 
= i 
0.4- Io fitting for total App, Arg(lAnl>1)  - 
[ and Arg(Mu>450 GeV) | 
0.3L J 


lo fitting for Apg(My>450 GeV) and Cy 


0-230 140 150 160 170 180 
Mz(GeV) 


FIG. 2. Parameter space scan of 1 o for all three asymmetries, Al, App(| An |> 1) and 


Arp(M > 450 GeV) is shown in red. The 1 o fitting for Arg(M;r > 450 GeV) and the total 


cross section Or. 


ea 5.7%-16.7% 27.5%-50.0% 23.1%- 74.5% |7.5 + 0.48 


TABLE I. 1 o benchmark points for all three asymmetry measurements. k-factor = 1.3, m}=172.5 


GeV for og. 


rapidity region while the selection cut are more efficient for the central events. Consequently, 
the cut efficiency at high invariant mass is quite low |20], which may further decrease the 
effective total cross section. Polarization of top quark in the events sample also effect the 


cut efficiency. 


2. W' 


To resolve the tension between cross section and total asymmetry in the Z’ model, the 
charged current process in t-channel may give better fit which has smaller interference effect 
due to the dd initial state. We plot the allowed parameter regions for the t-channel charged 
current model in in Fig Bl 
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FIG. 3. Parameter space scan of 1 ø for all three asymmetries, Al, App(| An |> 1) and 


Arg(M;r > 450 GeV) is shown in red. The 1 o fitting for Arg(M;r > 450 GeV) and the total 


cross section Or. 


The 1 o asymmetry region of all three measurements corresponds to a larger total cross 
section which is outside the 1 o fit of latest o,;r measurement. However, various efficiency 
effects discussed in the last paragraph of Z’ session may significantly reduce the measured 


cross section. 


3. Diquark 


We use the anti-triplet diquark that couples to tug to illustrate the feature. Similar to 
the W” case, there also exist diquark scalars whose couplings are of t°d@ and these diquark 


scalars contribute to dd — tt instead. 
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Figure [4] gives the 1 ø fitting for the anti-triplet diquark scalar with maximal flavor vio- 


lation. The 1 ø region also exists for the anti-triplet diquark scalar for all the measurements 
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FIG. 4. Parameter space scan of 1 ø for all three asymmetries, All, App(| An |> 1) and 


Arg(M;r > 450 GeV) is shown in red. The 1 o fitting for Arg(M;r > 450 GeV) and the total 


cross section Oyr. 


in asymmetries Afp. But the corresponding total cross section oj; are also larger than the 
measured value over 1 ø. In addition, the best-fit for cross section and the mass dependent 
asymmetry is over 1 TeV which makes the do/dM;; measurement very difficult to fit as 
shown in a The latest simulation by also showed that the tt events generated by 


diquark scalar had a higher cut efficiency at high energy therefore the anti-triplet diquark 
fitting is worse than the W”. 


II. IMPLICATIONS AT THE TEVATRON AND LHC 


After fitting the top forward backward asymmetries in different kinematical regions, we 
discuss the other Tevatron bounds for the models and the LHC predictions that can be soon 
tested in this section. 

The Large Hadron Collider (LHC) is a proton-proton collider with centre-of-mass energy 


7 TeV in the first two years running. Unlike at Tevatron where the axigluon effect only 


11 


appears as interference. The color octet axigluon of O(1 TeV) can be directly produced 
at the LHC and decay into dijet or tt. With significant decay branching ratio (BR) to tt, 
it provides additional handle to search it. The study of axigluon at the LHC has been 
performed by [7]. ATLAS collaboration has recently released the search for dijet resonance. 
The latest data has ruled out axigluon from 0.6-2.1 TeV by assuming axigluon coupling is 
only gs. The axigluon model in [5] has a even larger coupling comparing with the ATLAS 
paper and therefore, the model receive much more server constraint. 

For neutral gauge boson like Z’, the flavor violating vertex of ut will lead to large uu — tt 
or uu — tt scattering with Z’ exchange in the t/u-channel. The same-sign positive top quark 
pair (wu — tt) becomes particular interesting at the LHC given its large u-valence quark 
parton flux nt 

In addition, with large ut coupling, the tZ’ associate production is not negligible. Since 
Z' equally decays into ut and tu, the associated production tZ’ or tZ’ will contribute to 
tt+ j, tt+ 7 and tt+ j final states. Again, since the LHC is proton proton collider, the tt + j 
dominates the same-sign top production. The tt + 7 will appear in the inclusive tt search. 
Since the 1 o parameter space of all the asymmetry constraints corresponds to smaller tt 
pair production, the additional tt + 7 may in principle help to ease the tension at Tevatron. 
However, if it significantly contribute to the tt, the same amount of same-sign top quark 
will arise. 

Figure [5] (a) gives the pp — tt production rate at Tevatron and the 7 TeV LHC with the 
Z' in the 1 o fitting for all three asymmetry measurements. The pp — tt + tt at Tevatron 
is between 0.7-1 pb for these best fit points. CDF measured only 3 events for 2 fb~! 
with the acceptance range from 1.5% to 3%. The best fit points all predict 15-30 same-sign 
pure leptonic top events before selection cut but with one b-tagging. Even though these 
events from t-channel vector boson exchange may suffer from a low cut efficiency comparing 
to the t-channel light scalar exchange considered in Ref. [18], the Z’ model is strongly 
constrained by the same-sign top quark scattering. At Tevatron, the same-sign top due to 
tZ’ associate production is then much suppressed at the Tevatron due to significant phase 
space suppression. 

The uu — tt scattering get significantly enhanced at the proton-proton collider LHC. 
The production rate can reach 200 pb. Therefore, even at very early running of LHC with 


about 30 pb! data and requiring two b-tagging jet, the event number before kinematic cut 
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FIG. 5. (a) o(pp —> tt) at Tevatron and the 7 TeV LHC; (b) o(pp > tZ’+t+Z’) at Tevatron 
and the 7 TeV LHC. Both (a) and (b) are based on model parameter of Z’ in the 1 o fitting of all 


three asymmetry measurements as listed in Table [I] 


is about 70 and the same-sign top quark tt events is expected to be O(10). 
For W’ or diquark scalars with flavor violation, since they are electrically charged, it will 
only contribute to tt as at Tevatron. However, since the W’ or diquark & has a dui or 


ba 


With W’ and diquark scalars of typically above top quark threshold, they can decay into t 


ut coupling, the dg — tW’ or ug — td production is significant as shown in hd 


plus one hard jet. The signal is then tt plus one hard jet and should appear in the inclusive 
tt searches. The diquark case has already been calculated in our early paper [11]. 

Figure [6] gives the production of top quark plus W’ at hadron colliders. For My lighter 
than 400 GeV, the production rate is about 0.1-1 pb at the Tevatron. As we discussed 
in the previous section, the 1 o fitting parameter space for all the asymmetries constraint 
corresponds to the larger cross section region. |f| The new contribution to tt+ 7 will increase 
the tension between App and og. For heavy W’ above 400 GeV, due to large phase space 
suppression, the production rate at Tevatron can then be neglected. However, at the LHC, 
even with 7 TeV centre-of-mass energy, the production rate is O(10 pb). Even though 
Are that the single top production at Tevatron puts a strong constraint on the W’ — u —b coupling. 


However this constraint does not apply to general W’ models. 
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FIG. 6. pp > tW'>j+tt 


gg — tt dominates the tt at the LHC which does not interfere with W’, the tW’ itself may 


significantly increase the tt rate. 


II. INDIRECT CONSTRAINTS FOR MODELS 


A. Axigluon with flavor protection 


As shown in the previous section, only non-universal axigluon models can provide the 
positive asymmetry. Being a color octet with strong coupling strength, this GIM violation 


axigluon will then lead to significant flavor changing neutral current (FCNC) effect. 
L=igrgy" (HT); Prg TGS, + igray" Hr) Pap TG, (8) 


Flavor violation thus can arise from the non-universal gauge couplings due to the rotation 


between mass eigenstate and gauge eigenstate. 


UL = Vru;, dr = Vid, UR = VRur, dr = Vidr (9) 


The effective coupling in horizontal space is 


VARVE Vit eve, vel eve, vel HEVE, (10) 
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The rotation from mass eigenstate to gauge eigenstate for up and down type quark 
respectively is completely unmeasurable in the weak interaction. The only observable is the 
mixing in charge current transition which is categorized as CKM matrix i 

To avoid flavor violation in the down sector, one may introduce a U(1)a symmetry 
which acts only on the down sector with different eigenvalues for different generations, but 
does not distinguish the handedness of the quarks. Then the down-quark sector are diagonal 
with Ve = vg = Il so that there is no FCNC at all in the Bs, Ba or neutral K system. For 
simplicity, we take further the rotation matrix of right-handed up-quark sector as Vp = 1. 
Then one can explicitly determine the left-handed up-type quark rotation based on the 


known CKM matrix using veva = Voxm 
Vi = Vexm (11) 


Nowadays the Wolfenstein parametrization is widely used to express the CKM matrix 
in terms of four parameters (A, A, p and 7). To keep the unitarity of CKM matrix to all 
orders of A, we adopt in the following a definition of Wolfenstein parameters proposed in 


0l. Then the effective coupling between up quark and charm quark is 
(ViNHEVE J2 = —A?°(in — p +1). (12) 


Under the assumption of above rotations, the FCNC operators only arise in left-handed 


and mixing between first and second generation up-type quarks as 
1 =Q a) (a 1 Te 77 & 
a 0.34 Aanruch) + we Eu); (13) 


where the following decomposition satisfied by the color SU(3) fundamental representation 


has been implemented 


1 1 
Io). = dl — OapOne « (14) 
2 6 
Under Fierz transformation 
ra = (E eh), (15) 


the effective AC = 2 Hamiltonian can be expressed as 


HAG = C(u) (aye?) (Gi yuch) (16) 


T Tt was argued in i that the axigluon model in |] suffers serve bounds from Ba mixing. However, the 


calculation seems to be done by assuming both up and down quark sectors transform like CKM rotation. 
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and the leading order Wilson coefficient at the scale mg is 
2 44,10 “2 
g AtA (1 = p +i) 
C I. 17 
It means that the D? — D? mixing in this axigluon model has A!" suppression due to CKM 
rotation. The RGE running of the above Wilson coefficient is well known, 


ee 6/23 
C= ( — C(me) . (18) 


Notice that the hadronic matrix element of AC = 2 operator is 


n0|-a OL 2 
(D’jazy"ezaryuch|D?) = 3/ompBo(i) (19) 


Just like the B? — B° mixing case, one may define the renormalization group invariant 
parameter Bp by 

Bp = (a5(tte))~°/ Bote) , (20) 
which should be O(1). Then the axigluon induced AC = 2 effective operator contributes to 
mass difference of neutral D system as 


6/23 8r fAmMpBp as(me)A*d?((1 — p)? + n°) 


Amp = a;(mq) 
9 m3 


(21) 


The axigluon model can also induce AC = 1 effective operator which would in principle 


affect D? — D? mixing by 


M-T] ze iD mR eee, 5 
( 2 ), FA | eff | oe mp — En + ie (22) 
Actually the experimental observation of comparably large mass and width differences 
A AT 
o= P= 0.98t04%, ya 2 = (0.83 40.16)% (23) 
Tp 2p 


strongly implies that they are dominated by the long distance effects of the SM AC = 1 
operators. Therefore the axigluon induced AC = 1 terms could be safely neglected as they 
should be much smaller than the tree-level SM AC = 1 terms. 

Taking the Wolfenstein parameters as [32 


A=0.812, \=0.2254, p=0148, 71 =0.351 (24) 


8 The RGE running is actually dependent on Ny, the number of active flavor via do = 11 — 2Ny+/3. From 
the scale mg down to ue, Ny changes correspondingly from 6 to 4. But numerically this effect is small 


and we will simply take Ny = 5 in the RGE running. 
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and fp = 207 MeV , we obtain 


A {Tey \* 7B 
( ze) = 0.082% ( u ) a. (25) 
Pp axigluon Meg 1 


which is roughly one order of magnitude smaller than the experimental result. 


B. Electroweak constrains on the W’ model 


In general, the W’ must generate its mass through gauge symmetry breaking, then some 


other neutral component in the W’ symmetry breaking sector (for instance W} in the SU(2)r 


symmetry breaking) would inevitably mix with W? and some extra U(1) so that the W’~ 
would be charged under U(1)em. As a consequence, there is a large Z-Z’ mixing which is 
constrained by the electroweak precision test. In general, the bound from EWPT is subtle 
since different fermions W/Z boson couplings are modified in different ways which may even 
depends on models and a careful global fit is needed. The full results for the W’ model to 
explan the top forward backward asymmetry will be presented elsewhere. Since the overall 
modification for fermion Z boson coupling is small (except for some right-handed quarks 
charged under SU(2)r), we only consider the tree-level Z-Z’ mixing as a rough estimation. 
For observables that strongly depends on u/d/br-Z coupling, such as g, Qw(Cs), etc., their 
deviations from the SM results are still at the same level as Z mass which is transmitted 
into W boson mass. 

We can start to consider a simple SU(2)p x U(1)x x SU(2); model to estimate how 
large is the electroweak constraint for the Z’ — Z mixing. The SU(2)r is separated from 
SU(2), to avoid the troublesome W’ — W mixing. The two double Higgs field hz and hr 
are charged under SU(2); x U(1)x and U(1)x x SU(2)r respectively. The higgs fields 
get their vacuum expectation values (hy) = u, and (hr) = ur which spontaneously break 
SU(2)r x U(1)x x SU(2), into the diagonal group U(1)em. In order to raise the Z’ mass 
so we have less constrain from the Z-Z’ mixing, we choose hr transform as a triplet under 
SU(2)r so mg = V2myw:. The gauge quantum number for hz and hr are (0,1,1/2) and 
(1, 2,0) under SU(2)p x U(1)x x SU(2), respectively. For SM fermions, at least the quark 
doublet (t, d)r is charged under SU(2)r (It is possible to have some extra hidden fermions 
charged under SU(2)r and U(1)x to cancel the gauge anomaly). For the rest SM fermions, 


their quantum number is the same as the SM one if one replace their hyper charge as the 
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U(1)x charge. The quantum number for (t, d)r and (u, b)r under SU(2)pz x U(1)x x SU(2)L 
are (1/2, 1/3, 0). 
The kinetic term for the link fields Tr[(D,.hı)'(D,.hr)] + Tr[(D,.hr)'(D,hr)] becomes the 


mass terms for the massive gauge bosons. The mass matrix of the gauge bosons is 


r 29RUR —29RIXxUR 0 AR 
7 (az A At) —2grgxuh gluh + uz) —gxgiuz Ar l (26) 
0 —9x 9LUy, gru, A, 


We introduce the parameter € = u? /2u?, < 1 which shows that the right-handed symmetry 
breaking is only a perturbation. 
This matrix can be diagonalized by means of an orthogonal matrix which we shall call 


R: 


AR Ay 
AX | =R | Zz, |. (27) 
Ay Z, 


where the mass eigenstates are denoted by A, Z, and Z’. The eigenstate A is massless and 
identified as the photon. The couplings of our theory are related to the electric charge by 


e e e 


UX Sap? = 


28 
cos b cos Ow ` en 


IR = in dcos Oy ’ sin Ow 
where Oy is the weak mixing angle (in the limit € — 0) and d is an additional mixing angle. 


The other two eigenmasses are 


1 0 1 
2 2y 2 | 2 X ae 4 
= — =p 1 = = 1 Fa 5 29 
md, = Sub (ah +98) [1 | = zielen, 9 
my = Rokto [it erial = ot esinta], 60) 
2 (gR + 9x) 2 
where we have dropped O(e?) terms and 
1 1 1 
ae ee (31) 
I Gr Ix 


Clearly, Z is identified with the SM Z boson while Z’ is referred to as the heavy Z boson. 


For small e, the mixing matrix R has the following approximate form: 


sin ġ cos Ow cos œ cos Oy sin Ow 
33 2 ‚4 
R= | sindsindw + a cos ġ sin dw — g — cos Ow ; (32) 
— cos ġ + € cos ġ sinf ġ sin ġ + ecos? ọsin?’ ġ  —e cot Ow cos ọsin? & 
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2 
= 


0.0 ALAA AA ALVA | 
300 400 500 600 700 800 900 
Mz (GeV) 


FIG. 7. The excluded region on the parameter spaces of Wr model based on tree level contribution 
to T parameter alone from the new physics at 90% C. L. The green plus yellow region are is the 
excluded region for SM Higgs mass mp = 114GeV while the red plus yellow region is the excluded 
region for SM Higgs mass mp, = 150GeV. 


from which it is simple to derive the SM fermion couplings. The SM fermion and Higgs 


couplings to Z and Z’ can be written as 


grsin sin wT? + gx cos dsin dwQx — gr cos OwT? 
= ——— (sin? OwQ — T? + esin? ġ cos? oT} — esin* oQx) : (33) 
sin Oy cos Oy 
gr(—cos@ + ecos dsin? o)T} + gx(sind + ecos? osin? 6)Qx — gr(-e cot Ow cos dsin® ¢)T? 


. . 3 
= e E ae EN AE eg Oy) i (34) 
cos Ow sin cos @ sin“ Oy 


In the limit of large SU(2)r breaking vev (e < 1) and small mixings (® — 0), the Higgs 


current can be approximated as (we drop the € cos? ¢sin® & term) 


1 
Joh) = —59x sin b(h'DYh) + hc. , (35) 
which induce a dimension six operator 


2 2 
-IRT IK sint (nt Dh)? . (36) 


anOn = 2m 
ZG 
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which coincident with Eq. that Am?, = —esin ¢*m?. We can calculate the correspond- 
ing T parameter from the tree level gauge boson mixing, 


2 ‚4 
aRU esin @ 
T = — = 


a 7 (37) 


Using the SM model mz, Gp (the life time of T) and a = e?/4r as the basic input param- 
eter, we can calculate the allowed parameter space including the Higgs radiative corrections 
according to the most recent results: S = 0.03 + 0.09 and T = 0.07 + 0.08 (with 87% strong 
correlation) hal, The results are presented in Fig[7] We can see that for sufficient heavy Z’ 
and strong coupling gr (for instance, gr = 2, mz = 900 GeV which is used in Ref. lod), it 


is well above the excluded region. 


IV. CONCLUSIONS 


We discuss the feature of various models for the top quark forward-backward asymmetry 
anomaly at Tevatron, using the latest CDF measurements on total asymmetry in lab frame 
A, the rapidity dependent asymmetry Arg(| An |> 1), the mass depedent asymmetry 
Arg(M;r > 450 GeV) and the total tt production cross section oy. 

The axigluon model in [5] has difficulty in explain the large rapidity dependent asymmetr 
and the mass dependent asymmetry simultaneously. In addition, the latest dijet search 
at ATLAS has ruled out the parameter region that is relevant to top Arg. On the other 
hand, in contrast to the conclusion in Ref. sl, the model itself does not suffer from the 
flavor Bg mixing under flavor protection U(1)q and a careful calculation shows that their to 
D? — D? is still one order lower than the current experimental bound. 

The t-channel Z’ [9], W” and anti-triplet diquark models all have parameter re- 
gions that satisfy all three asymmetry measurements within 1o. However, the corresponding 
production cross section predicted by the 1 o asymmetry requirement in the Z’ model are 
always significantly below the 1 ø of cross section measurement. The best fit point of Z’ 
is about 700 GeV with purely righthanded coupling g > 1.8 which corresponds to 6.9 pb. 
However, this best fit point will generate a large number of same-sign top quark events at 
Tevatron which is at least five times larger than the SM prediction. We conclude that the 
Z' model is very difficult to be consistent with all the measurements. 

Both W’ and anti-triplet diquark models predict the cross sections are larger than the 


measurement but various factors can lower the survival efficiency after cuts in these models 
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to ease the tension between asymmetry and cross section. The best fit point for anti-triplet 
diquark lies in very high mass region and with better survival efficiency Id. it is difficult to 
fit the differential cross section do/dM,. A rough estimation for W’ model shows that the 
bounds from electroweak precision tests are weak due the heavy Z’ and strongly coupling 
gr. Therefore, we conclude that the best model is the t-channel W’ model at the current 
stage. To test such model directly, we also use the 1 o asymmetry parameters to compute 
the production rate of tt + j from tW” at 7 TeV LHC and the production rate is 10%-50% 
of SM tt. 

Last, we want to mention that the latest NNLL calculation olm, = 173.1 GeV) = 
6.30 + 0.194233 pb [1] is significantly lower than the experimental results. If the result does 
not significantly change for m; = 172.5 GeV which is used for Tevatron experiments, then 
the fits for t-channel W’ and anti-triplet would be better while the t-channel Z’ would be 


worse. 
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NOTES ADDED 


While this work was being delayed by the huge earthquake in Japan, Ref. appeared, 


which overlap with ours in the study of fitting different models on different top forward 
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backward asymmetries. Our results agree with quantitatively with theirs for fitting different 
top forward backward asymmetries and the total tt cross section. However, we notice that 
the axigluon model and the heavy Z’ model are severely constrained by the dijet search at 
the LHC (ATLAS) and the same sign dilepton search at the Tevatron (CDF). Therefore, we 
conclude that the heavy W’ model is the most promising one at present. We also consider 


the indirect bounds for different models. 
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